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Abell48 - a rare WN-type central star of a planetary 
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ABSTRACT 

A considerable fraction of the central stars of planetary nebulae (CSPNe) are hydrogen- 
deficient. Almost all of these H-deficient central stars (CSs) display spectra with strong 
carbon and helium lines. Most of them exhibit emission line spectra resembling those 
of massive WC stars. Therefore these stars are classed as CSPNe of spectral type [WC] . 
Recently, quantitative spectral analysis of two emission-line CSs, PB 8 and IC 4663, re- 
vealed that these stars do not belong to the [WC] class. Instead PB 8 has been classified 
as [WN/WC] type and IC4663 as [WN] type. In this work we report the spectroscopic 
identification of another rare [WN] star, the CS of Abell48. We performed a spectral 
analysis of Abell48 with the Potsdam Wolf-Rayet (PoWR) models for expanding at- 
mospheres. We find that the expanding atmosphere of Abell48 is mainly composed of 
helium (85 per cent by mass), hydrogen (10 per cent), and nitrogen (5 per cent). The 
residual hydrogen and the enhanced nitrogen abundance make this object different 
from the other [WN] star IC 4663. We discuss the possible origin of this atmospheric 
composition. 

Key words: stars: abundances stars: AGB and post-AGB stars: mass-loss 
stars: Wolf Rayet planetary nebulae: general planetary nebulae: individual: 
PN G029.0+00.4 



1 INTRODUCTION 

During their evolution, low mass stars in a certain mass 
range become central stars of planetary nebula? (CSPNe). 
Central stars (CSs) are post-AGB stars with effective tem- 
peratures above 25 kK, hot enough to ionize circum-stellar 
material, which was shed off during the AGB phase. The 
ionized circum-stellar gas becomes visible in the optical as 



* Based on observations obtained with the South African 
Large Telescope (SALT), commissioning programme 
2010-1-RSA_OTH-001 and programme 2011-3-RSA_OTH-002. 
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a planetary nebula (PN) around the CS. While most of the 
low mass stars stay hydrogen-rich at the surface through- 
out their life, a fraction of about 20 per cent of CSs show a 
hydrogen-deficient surface composition. Almost all of these 
H-deficient CSs exhibit spectra with strong carbon and he- 
lium lines, indicating that their surface is composed pre- 
dominantly by these elements (e.g. Gorny fc Tylend"a||2000| 
Werner & Heber|1991l. 



About half of the H-deficient CSs show emission line 
spectra resembling those of massive WC stars, i.e. carbon- 
rich Wolf-Rayet (WR) stars with strong stellar winds. In 
analogy, the CSs with bright carbon emission lines are clas- 
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sified as [WC] stars, where the brackets should distinguish 
them from their massive counterparts. 

Recently, quantitative spectral analysis of two CSs with 
broad emission lines revealed that they do not belong to 
the [WC] class. The first one was PB 8 ( |Todt et al.||2010a[ ), 
which should be classified as [WN/WC] type, and the other 
one was IC4663 ( |Miszalski et al.|2012[ ), which shows a [WN] 
type spectrum. In both cases the atmospheres of these stars 
are mainly composed of helium. The atmosphere of PB 8 
also contains hydrogen and traces of carbon, nitrogen, and 
oxygen, while IC 4886 has an almost pure helium atmosphere 
with only traces of nitrogen. 

For some other CSs the membership to the [WN] class 
was discussed, but could not be proven | |Todt et al.|2010bp . 
DePew et aL] ( |2011| mentioned the CS of Abell48, which 
shows a WN-type spectrum, but it was not clear whether 
this is really a CSPN and not a massive WN star with a 
ring nebula. In this paper, we present a quantitative anal- 
ysis of the optical spectrum of the CS of Abell48, which 
reveals indeed a WN-like atmospheric composition. We will 
also present arguments for the low-mass nature of Abell48. 

The paper is organized as follows: Section [2] gives an 
overview of previous observations of Abell 48. In Section [3] 
we describe our observations and the data reduction. A short 



Table 1. The central star of Abell 48. 



4] The spectral 
deals with the 



description of our methods is given in Section 
analysis is presented in Section [5] Section [6 
nebula, while the results and implications are discussed in 
Section [3 



2 ABELL 48 



Abell 48 was discovered by Abell ( 1955 1 on the original 



plates of the Palomar Sky Survey. In the red photograph 
from this survey (presented in fig. 1 in |Abell|19 66), Abell 48 
appears as a slightly elliptical double-ringed nebula with a 
major axis diameter of about 40" (see also the left panel 
of Fig. [T] for the Digit ized Sky Survey II (DSS -II) red band 



( McLean et al.||2000[ ) image of the nebula). |Abell| ( |1955 
19661 classed the nebula as PN (listed in the discovery pa- 
pers under numbers 36 and 48, and named in the SIMBAD 
database PN A55 36 and PN A66 48, respectively), while 
its CS has not been classified. Since that time, the PN clas- 
sification of Abell 48 was generally accepted and the nebula 
was included in the Strasbourg-ESO Catalogue of Galac- 



1992 1 under the name 



tic Planetary Nebulae (Acker et al. 
PN G029.0+00.4. 

Abell 48 is a strong radio source. With the 1.4 GHz flux 
of 159±15mJy (Condon & Kaplan 19981, it is one of the 



brightest objects among the PNe covered by the NRAO 
VLA Sky Survey (NVSS; |Condon et al.|p98l). The circu- 



lar shell of Abell 48 was clearly resolved in the Multi- Array 



Galac tic Plane Imaging Survey (MAGPIS; Helfand et al. 
2006), carried out with the Very Large Array (VLA)Q In 
the MAGPIS 20 cm image Abell 48 appears as a thick shell 
of the same size as the optical nebula (see Fig. [T]). 

Abell 48 was also detected in several major mid- and 
far-infrared (IR) surveys, while its ring-like shell has been 
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resolved in the IR only with the advent of the Spitzer Space 



Telescope (Phillips & Ramos-Larios 2008 Wachter et al 



20101. Figure [1] includes the Spitzer 8 and 24 fim images of 



Abell 48 obtained with the Infrared Array Camera (IRAC; 
|Fazio et al.||2004| ) and the Multiband Imaging Photometer 
for Spitzer (MIPS; |Rieke et al~1|2004| | within the framework 
of the Galactic Legacy Infrared Mid-Plane Survey Extraor- 
dinaire (GLIMPSE; |Benjamin et~ldT1|2003"| ) and the 24 and 
70 Micron Survey of the Inner Galactic Disk with MIPS 
(MIPSGAL; |Carey et aL]|2009[ ), respectively. In the 8 fim 
image the nebula shows a bright inner shell and a diffuse 
outer halo with an angular extent equal to that of the outer 
optical shell, while at 24fj,m its morphology is very similar 
to that of the radio nebula. 

To our knowledge, there was not any significant study 
of Abell 48 until recently. The interest in Abell 48 has arisen 



after Wachter et al. (20101 carried out optical spectroscopy 
of its CS and classified it as WN6. Based on this result, 



DePew et al. (20111 and Bojicic et al. (20121 have specu- 



lated that Abell 48 might be another member of the [WN] 
or [WN/WC] class (|Todt et aL]|2010a[), but there has not 



been a detailed study yet. 

Table[l]gives the coordinates and photometry of the CS 
of Abell 48. The B and R magnitudes are from the Guide 
Star Catalogue (GSC2.2J5 the / magnitude is from the 
DENIS (Deep Near Infrared Survey of the Southern Sky) 
database ( |DENIS Consortium]|2005| >, the J, H, K s magni- 
tudes are from the 2MASS (T wo Micron All Sky Survey) All- 
Sky Catalog of Point Sources ( |Cutri et al.|2003 1, and the co- 
ordinates and the IRAC magnitudes are from the GLIMPSE 
Source Catalog (I + II + 3D^] 



3 OBSERVATIONS AND DATA REDUCTION 

Abell 48 was observed within the framework of our ongo- 
ing programme of spectroscopic follow-up observations ded- 
icated to evolved massive star candidates that were re- 
vealed through the detection of their IR circumstellar nebu- 



lae ( Gvaramadze et al. 


2009 2010a|b|c 2011 2012 Stringfel- 


low et al.||2011 2012 


Burgemeister et al. 20131. Although 


Abell 48 was not included in our (Gvaramadze et al.|2010al 



1 The authors of this survey classed Abell 48 as a high-probability 
supernova remnant candidate. 



http: / / vizier.u-strasbg.fr / viz-bin /VizieR?-source=I/271 



http://cdsarc.u-strasbg.fr/viz- bin/CatYli/293 
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Figure 1. From left to right: DSS-II red band, Spitzer 8 and 24 /im, and the VLA 20 cm images of Abell 48 and its central star (marked 
by a circle in the DSS-II and 8 fim images). The orientation and the scale of the images are the same. 



list of compact nebulae with central stars (because it was 
identified in the SIMBAD database as a PN), we decided to 
observe it to check whether its CS is an ordinary WN star 



(as suggested by Wachter et al. 2010 1 or we deal with an 
extremely rare object - a [WN]-type CSPN. 

Spectral observations of Abell 48 were obtained with the 



Southern African Large Telescope (SALT; Buckley et al 



2006 O'Donoghue fc et al. 2006) on 2011 August 25 dur- 



ing the Performance Verification phase of the Robert Stobie 



Spectrograph (RSS; Burgh et al. 2003 Kobulnicky et al 



2003), and one more spectrum was taken on 2012 April 12. 



The long-slit spectroscopy mode of the RSS was used, with 
a slit width of l'.'25 and a position angle of 0° for all ob- 
servations. The seeing was stable in the range of l'.'6 to 
2'.'0. We utilized a binning factor of 2, to give a final spatial 
sampling of 0'.'258 pixel" 1 . The Volume Phase Holographic 
(VPH) grating GR900 was used in these observations to 
cover a total spectral range of 4300-7300 A with a final re- 
ciprocal dispersion of ~0.97A pixel -1 and a spectral reso- 
lution FWHM of 5-6 A. Three exposures were taken, each 
with 15 minutes. Spectra of the Xe comparison arc were ob- 
tained to calibrate the wavelength scale. Spectrophotometric 
standard stars were observed during twilight. 

RSS has three mosaiced 2048x4096 CCDs as a detec- 
tor. First, data for each CCD were overscan subtracted, 
trimmed, gain and cross-talk corrected, and finally mo- 
saiced. This primary data reduction was done with the SALT 
science pipeline ( |Crawford et ai"||2010| ). Second, primary 
corrected and mosaiced long-slit data were reduced in the 



way described in Kniazev et al. (2008b I. Finally, all two- 



dimensional spectra were averaged. The spectrum of the sur- 
rounding nebula was subtracted using the IRAF background 
task and a one-dimensional spectrum of the CS was ex- 
tracted with the apall task. The two-dimensional spectrum 
of the surrounding nebula was averaged along the slit. 

SALT is a telescope with a variable pupil, and its illumi- 
nating beam changes continuously during the observation. 
This makes an absolute flux calibration impossible, even us- 
ing spectrophotometric standard stars or photometric stan- 
dards. Therefore the spectrophotometric standard stars were 
used only for a relative flux calibration, and the absolute flux 
level was scaled to the photometric R magnitude. 



4 METHODS 

4.1 Spectral modelling 

For the analysis of the spectrum we used the Potsdam 
Wolf-Rayet (PoWR) models for expanding atmospheres (see 
|Grafener et al.|2002|[Hamann fc Grafener|2004[ ). The PoWR 
code solves the non-LTE radiative transfer problem in a 
spherically expanding atmosphere simultaneously with the 
statistical equilibrium equations while accounting for energy 
conservation. Iron-group line blanketing is treated by means 
of the superlevel approach. Wind inhomogeneities in first- 
order approximation are taken into account, assuming small- 
scale clumps. In this work, we do not calculate hydrody- 
namically self-consistent models, as in |Grafener fc Hamann| 
1 2005 ) , but assume that the velocity field follows a /3-law 
with /3 — 1 with the mass-loss rate as a free parameter. Our 
computations include complex atomic models for hydrogen, 
helium, carbon, nitrogen, oxygen, and the iron-group ele- 
ments. Within our group, this code has been extensively 
used for quantitative spectral analyses of massive WR stars 
(e.g. Grafener et al.||2002 Liermann et al. |2010| ) and WR- 



type CSPN (e.g. |Todt et ai.||20i"0a I 



4.2 Spectral fitting 

The typical emission-line spectra of Wolf-Rayet stars are 
predominantly formed by recombination processes in their 
dense stellar winds. Therefore the continuum-normalized 
spectrum shows a useful scale-invariance: for a given stellar 
temperature T* and chemical composition, the equivalent 
widths of the emission lines depend only on the ratio be- 
tween the volume emission measure of the wind and the area 
of the stellar surface, to a first approximation. An equivalent 
quantity, introduced by |Schmutz7 Hamann, & Wcsso lowski| 
( |1989[ ), is the transformed radius 

-I 2/3 



Rt — R* 



2500 kms- 



10 



' 4 M yr- 



(1) 



Different combinations of stellar radii R, and mass-loss rates 
M can thus lead to nearly the same emission-line strengths. 
In the form given here, the invariance also includes the 
micro-clumping parameter D, which is defined as the den- 
sity contrast between wind clumps and a smooth wind of the 
same mass-loss rate. Consequently, mass-loss rates derived 
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empirically from fitting the emission-line spectrum depend 
on the adopted value of D. The latter can be constrained by 
fitting the extended electron scattering wings of strong emis- 



sion lines (e.g. Hillier f991 Hamann fc Koesterke|1998 1 



5 ANALYSIS 



5.1 Stellar parameters 



In order to find the model that best matches to the obser- 
vation, we choose a systematic approach: We begin with a 
first determination of log T„ and log R t with the help of the 
already existing mo del grid for massive WN sta rs of the late 
subtypes (WNLQ jHamann fc Grafener 20041, which was 
recently updated with improved atomic data. 

For each of the models in this grid, we calculated the 
reduced xt 



1 N 
\ 2 = — V 



f? od ) 2 



(2) 



f° B denotes the by-eye rectified observed spectrum and 
jmod ^ e con tinuum normalized synthetic spectrum of the 
model, each with N data points. We consider R t and T* as 
the fitting parameters. We adjusted the uncertainty a so, 
that the reduced xt of the best fitting model equals unity, 
corresponding to a mean deviation of about 20 per cent. 
The best fitting model is the one with logT»/K = 4.85 and 
logRt/Re = 0.9. 

The models of this grid had been calculated with a 
bolometric luminosity of logL/L0 = 5.3 and a density con- 
trast of D = 4, as appropriate for many massive WN stars 
(Hamann et al. 2006[ ). However, in SectionJTJwe present argu- 
ments that Abell48 is in fact a PN. Therefore we calculated 
a grid of models with parameters that are more appropri- 
ate for CSPNe, i.e. for the stellar luminosity and mass of 
L = 6000 L e and M = 0.6 M (see e.g . 
20051 iMiller Bertolami & Althaus|2007|) 



Schonberner et al. 



Moreover, with our first estimate of R t and T» we calcu- 
lated a series of models with different chemical abundances 
to determine the most appropriate atmospheric composition. 
The best fit is achieved with the abundances presented in 
Table [3] as described in Section [5.2| These abundances were 
used for our [WN] grid. 

Since there is no UV observation available that would 
contain spectral lines with P Cygni profiles, we had to 
deduce the wind terminal velocity, v x from the width of 
the optical emission lines. The best fit is achieved with 
Voo = f000 ± 200 km s -1 . Line broadening by microturbu- 
lence is also included in our models. From the shape of the 
line profiles we estimated the microturbulence to be about 
100 kms -1 . 

We repeated our y^-Ht analysis for this grid of [WN] 
models to get the uncertainties in _R t and T* (see Fig. [2J. 
The lcr-uncertainties can be obtained from the contour with 
xt = 2, i.e. where the value of xt nas increased by 1 com- 
pared to the best-fitting model, for which xl — L The best- 
fitting model from this [WN] grid, has T* = Tll 1 ^ kK and 
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Figure 2. Contours of the reduced \t f° r our grid of [WN] mod- 
els. The open circles indicate the calculated models. Between 
these data points the contour lines are interpolated. Note that 
a = 1 corresponds to a \a confidence interval. Thus, the 

contour with \i — 2 indicates the la uncertainty of the parame- 
ters Rt and T* . Our final model is indicated by the red square. 
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Figure 3. Contours of the ratio between the line peak heights of 
Niv 7100 / Nin 4643 (green dashed lines), Niv 7100 / Nv 4933 
(blue filled ellipse), and Hell 5412 / Hel 5876 (black solid lines). 
The thick contours represent the measured values, the thinner 
lines the ltr uncertainty. The other symbols are the same as in 

Fig-0 
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Table 2. Ratios between the peak heights: measured ratios, in- 
ferred quantities, and ratios from the final model. 



Line ratio 


Observed 


Final model 


HenA5412 / He I A5876 


2.45±0.25 


2.2 


NivA7120 / NmA4634 


3.3±0.5 


2.3 


NivA7120 / NvAA4933/43 


2.8±0.3 


2.1 



\ogRt/RQ = 0.9 ± 0.2, where the given error margin corre- 
sponds to a doubling of the deviation between model and ob- 
servation, as measured by the reduced xl ( see Equation [2]). 

For a further refinement, we calculated the line-strength 
ratios between significant lines of He 11 / He I, N IV / N ill, and 
N IV /N v for this grid and compared them to the observation 
(cf. Table [2} . We used line ratios instead of absolute line 
strengths to diminish the influence of chemical abundances. 
In practice, it is often difficult to measure the equivalent 
widths of lines correctly, e.g. due to line blends. Therefore 
we used the peak height as a measure of the line strength. 

From our optical observation we estimate the uncer- 
tainty in the normalized continuum to be of the order of 10 
per cent in the region of the N 111 A4643 line and 5 per cent 
at the other diagnostic lines. We consider this as the uncer- 
tainty in the peak measurement, which we then use to infer 
a 10 per cent uncertainty in the measured He 11 / He I and 
N IV /N v peak height ratio and a 15 per cent uncertainty in 
the measured Niv /Nin peak height ratio. 

In Figure [3] we show contours of the N IV / N ill and 
He 11 /He 1 line peak ratios. The bold lines correspond to 
the observed values. These contours have two intersections. 
As a further criterion we now consider the N IV /N v ratio. 
Unfortunately, all models of the grid give a smaller value 
than observed, but the two models encircled by the shaded 
area in the upper left corner of the diagram come closest. 
After manual inspection of the spectral fits we finally de- 
cided for the final model with logRt/Ro = 0.851q'15 and 
T„ = 70i* kK, indicated by the red square in Figure pi The 
given uncertainties were estimated from the la contours in 
this figure. 

Remarkably, the final parameters derived under as- 
sumption of a CS luminosity hardly differ from those ob- 
tained above from the massive star grid. This demonstrates 
that the scale-invariance for models with same R t holds over 
a wide parameter range. 

Together with the normalized spectrum we obtained 
the synthetic spectral energy distribution (SED) in abso- 
lute units from the PoWR models. The synthetic SED was 
fitted to the flux-calibrated spectrum and photometric mea- 
surements by adjusting the distance and the reddening pa- 
rameter Eb-v (cf. Fig. [5|. The photometric magnitudes for 
the Johnson and 2MASS measurements were converted to 
fluxes with help of the calibrations by Holbcrg & Be rgeron] 
(20061, the IRAC measurements with help of the calibra- 



tions by Cohen et al. (20031. 

The best SED-fit was obtained with a color excess of 
Eb-v = 2.10 mag and the standard total-to-selective ab- 
sorption ratio Rv = 3.1 ( Fit zpatrick| 1999] ) . We estimated a 
distance of d = 1.9 kpc towards Abcll48 with the adopted 
stellar luminosity of 6000 Lq. 



5.2 Element abundances 

In an iterative process with the determination of T* and 
Rt, we varied the abundances of He, H, C, N, and O. As 
there is no UV-spectrum of the CS of Abell 48 available that 
would help to determine the abundances of the iron-group 
elements, we kept latter fixed to solar values. 

Helium. From the strengths of the helium emission lines, 
especially Hen 4-3, it is already obvious that the wind is 
mostly composed of helium. 

Hydrogen. None of our models is able to reproduce the 
Ha/He 11 and H/3/Hell blends together with the unblended 
Hell 7-4 line simultaneously as observed (cf. Fig.[6|. H/? fits 
best with a hydrogen mass fraction of 20 per cent, while 
the Ha/He 11 blend matches the observation with about 3 
per cent of hydrogen. As a compromise we adopt a value of 
10 per cent of hydrogen. A better resolved spectrum with 
a higher signal-to-noise ratio (S/N) would help to pin down 
the exact hydrogen abundance. 

Nitrogen. The nitrogen lines in the spectrum of the CS 
of Abell 48 are much stronger than in WN spectra of stars 
with similar parameters. We had to increase the nitrogen 
abundance to 5 per cent by mass to match the observed 
line strength of the Nv AA4933/4944 multiplet and most 
of the other nitrogen lines. The Niv A7125 line and the 
N iv/N v blend at 4600 A are unfortunately not well matched 
by our model. Models with only An = 3 per cent provide 
a much better fit to the NiiiA4634 line, while the strength 
of the Niv A7125 line is best reproduced by models with 
An = 7 per cent. Generally, it is hard to find a sufficient 
fit for spectral lines of three different ions from the same 
element. This might be an indication for an inhomogeneous 
wind. 

Carbon. In the absence of other carbon lines within 
the range of our observed spectrum we have to rely on the 
C iv A5800 line only. This line can be best fitted with a 
model that has a carbon abundance of 0.3 per cent by mass. 
The given uncertainty of ±0.1 per cent accounts only for the 
uncertainty of the continuum normalization 

Oxygen. There are no strong oxygen lines in the spec- 
tral range of our observation. A model with solar oxygen 
abundance of Xo — 0.6 per cent by mass predicts a weak 
emission line from O ill A5592 with a line strength that is al- 
most below the detection limit. A higher oxygen abundance 
would not be compatible with the observation. 

The final line fit is shown in Fig. [4] and the model pa- 
rameters are compiled in Table [3] 



6 THE NEBULA OF ABELL 48 

6.1 Position— velocity diagram along the slit 

In Fig. [7] we present the Ha intensity and radial velocity 
distributions along the slit, and the DSS-II red band image 
of Abell 48 for comparison. The extension of the nebula in 
the Ha line is very similar to that in the DSS-II image. Ha 
emission appears everywhere within the nebula and its peaks 
are clearly correlated with the CS and the shell's position. 
We measured the diameter of the nebula of 43" at a ten 
per cent level of the peak intensity in the Ha line. At the 
distance of Abell 48 of 1.9kpc (see Sections 5.1 and 7.3l the 
radius of the nebula is about 0.2 pc. 
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Figure 4. Optical spectrum of the CS of Abell 48: Observation (blue solid) vs. synthetic spectrum (red dashed). The observed spectrum 
has some gaps and gets much noisier towards shorter wavelengths. 




4.2 4.4 
log X [A] 

Figure 5. Observed flux distribution of Abell 48 (blue/noisy) in absolute units, including the calibrated spectrum and the photometric 
measurements (open blue squares, see TableJTJ, compared to the emergent flux of the model (red/smooth line). The model flux has been 
reddened and scaled to the distance according to the parameters given in Table [3] 



We calculated the line-of-sight velocity distribution 
(vhei) along the slit using the method and programs de- 



scribed in Zasov et al. (20001. To exclude systematic shifts 
originating from the known RSS flexure, the closest bright 
night sky line [0 1] A6363 A was used. Finally, only those 



velocity estimates where used which satisfy the criteria 
S/N > 3 and a v < lOkms" 1 . 

Fig. [7] shows that Vhe\ along the slit spans a rage of ~ 
to 50kms _1 . This range is very similar to that found for PNe 
(see e.g., |Torres-Peimbert et al.|2009 l. We measured t>hci = 
50.4 ± 4.2 km s~ in the area which is closest to the position 
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Figure 6. Test for the hydrogen abundance in the CS of Abell48. The observed profiles (blue/thick solid line) are compared with 
models (red/thin lines) for different hydrogen abundance, but otherwise the same parameters as of our best-fit model (see Table [3}. The 
tested hydrogen mass fractions are zero (dotted lines), fO per cent (thin solid lines) and 20 per cent (dashed lines), respectively. 



Table 3. Abell 48: Stellar and wind parameters. 



T, 


70 ±t 


kK 




1000 ± 200 


kms 


logflt 


o.85t :i 5 3 


i? Q 


logM 


-6.4 ±0.2 


MQyr" 1 


R* 


0.54 ±0.07 


Rq 


D 


4 


(density contrast) 


Eb-v 


2.10 ±0.05 


mag 


""i-ad 


50 ±4 


kms -1 


H 


10 ± 10 


per cent by mass 


He 


85 ± 10 


per cent by mass 


C 


0.3 ± 0.1 


per cent by mass 


N 


5.0 ± 2.0 


per cent by mass 


O 


< 0.6 


per cent by mass 


Fe 


1.4 x 10" 3 


per cent by mass 



Table 4. Line intensities of the nebula Abell 48. 

A (A) Ion F(A)/F(H/3) I(A)/I(H/3) 



of the CS and suggest that this velocity corresponds to the 
velocity of the mass center. With our spectral resolution 
(FWHM = 5.5 ± 0.5 A) for the range of the Ha line, we are 
not able to measure the expansion velocity directly. 



6.2 Physical conditions and chemical abundances 

The spectrum of the nebula of Abell 48 has never been 
studied before, because it is relatively faint and highly red- 
dened. An objective prism spectrum was taken by|Sanduleak 



DePew et al. 



(19761, where only the Ha line was detected 
( 2011 1 studied the CS of Abell 48 with the Siding Spring 
Observatory 2.3m telescope, but they did not publish any 
spectrum. Even with SALT we were not able to detect any 
useful signal in the spectral range below 4300 A. Our final 
ID spectrum of the nebula is shown in Fig. [8] 

Emission lines in the spectrum of the surrounding neb- 
ula were measured applying the MIDAS programs described 
in detail in |Kniazev et al.| ( |2004| |2005[ ). Table [2] lists the 
measured relative intensities of all detected emission lines 
relative to H/3 (F(A)/F(H/3)), the ratios corrected for the 
extinction (I(A)/I(H/3)), and the derived extinction coeffi- 
cient C(H/3). The latter corresponds to Eb-v of 2.15 mag, 
which is very similar to the value found independently from 
the spectrum of the CS (see Table [3|. 

The spectrum of the nebula was interpreted by the tech- 



4340 


H7 





170±0 


009 


0.451±0 


024 


4363 


[0 in] 





014±0 


005 


0.034±0 


013 
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1 


000±0 
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1.000±0 
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4959 


[O in] 


1 
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1.005±0 
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[0 in] 


4 
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018±0 


002 


0.004±0 
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He 1 





974±0 


024 


0.206±0 
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[O 1] 





042±0 


003 


0.005±0 
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6548 
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3 


103±0 
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0.282±0 


012 


6563 


Ha 


32 


555±0 
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2.906±0 
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6584 


[N 11] 


8 


831±0 
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0.770±0 
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6678 


He 1 





616±0 


015 
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761±0 


018 
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002 
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[S 11] 
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0.068±0 


002 


7065 


He 1 





434±0 


012 


0.022±0 


001 


7136 


[Ar in] 


2 


117±0 


051 


0.102±0 


003 


7237 


[Ar iv] 





057±0 


003 


0.003±0 


000 


7281 


He 1 





170±0 


010 


0.007±0 


000 


7320 


[O 11] 





174±0 


014 


0.007±0 


001 


7330 


[O 11] 





156±0 


014 


0.006±0 


001 


C(H/3) dcx 






3.16±0.03 




E B - 


V mag 






2.15±0.01 





nique of plasma diagnostics in the way described in detail in 



Kniazev et al. (2008b I. The electron temperatures T e (Olll), 
T e (N 11), were calculated directly using the weak auroral lines 
of oxygen [Om] A4363 and nitrogen [Nil] A5755. Both tem- 
peratures and the number density n e (Sn) are given in Ta- 
ble [5] From the detected emission lines we were able to de- 
termine the total element abundances for O, N, Ar, and He. 
These results are shown in Table[5]together with solar abun- 
dances ( jAsplund et al.|2009| ). 

The derived abundances of the a-elements O and Ar 
are lower than solar by about — 0.3dex. In contrast, He and 
N show abundances close to the solar values. 

The IR fluxes at wavelengths longer than 8 /im, which 
have been measured by Midcourse Space Experiment (MSX) 
satellite ( Price et al.|2001 1 and the Wide- field Infrared Sur- 
vey Explorer ( WISE; |Wright et al.|2010 l, show a strong ex- 
cess compared to the predicted stellar continuum. This ex- 
cess could be due to the low resolution of these instruments, 
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Table 5. Elemental abundances in the nebula of Abell48. 



Figure 7. Ha intensity and radial velocity distribution along 
the slit. Upper panel: DSS-II red band image of Abell48 with the 
position of the 1''25 slit shown by dashed lines. Middle panel: The 
relative flux of the Ha-line along the slit after continuum subtrac- 
tion. N-S direction of the slit is shown. The vertical line at r = 
corresponds to the position of the CS. The region rrfc 2" is shown 
with vertical, dashed lines to mark the area, where the average ve- 
locity i>hel = 50.4±4.2kms —1 was calculated. Bottom panel: The 
velocity profile of the Ha line, corrected by using the night-sky 
line [O i] A 6363 A. The horizontal line indicates the calculated ve- 
locity 50.4kms — 1 , that has to be close to the heliocentric velocity 
of the CS. 



which results in confusion between the stellar source and the 
nebular emission. 



7 DISCUSSION 

As discussed in the introduction, there were doubts whether 
Abell 48 is a PN with a low-mass central star or a ring nebula 
around a massive star. Now we discuss several arguments in 
favour of the PN status of this object. 



Quantity 



Abell 48 



Sun a 



T e (OIII)/K 

T e (NII)/K 
rie(SII) / cm" 3 

O+/H+(xl0 5 ) 
O++/H+(xl0 5 ) 
O/H(xl0 5 ) 
12+log(0/H) 

N+/H+(xl0 7 ) 

ICF(N) 

N/H(xl0 5 ) 

12+log(N/H) 

log(N/0) 

Ar++/H+(xl0 7 ) 

ICF(Ar) 

Ar/H(xl0 7 ) 

12+log(Ar/H) 

log(Ar/0) 

He/H 

12+log(He/H) 



11870±1640 
7200±750 
1000±130 

19.97±5.89 
6.56±2.74 

26.53±6.49 
8.42±0.11 

416.9±151.6 
1.55 
6.48±2.36 
7.81±0.16 
-0.61±0.19 

6.06±2.53 

1.21 
7.33±3.05 
5.86±0.18 
-2.56±0.21 

0.133±0.005 
11.12±0.02 



8.69 



7.83 
-0.86 



6.40 
-2.29 



10.93 



a Solar abundances from [Asplund et al.| |j2009) 

7.1 Abell 48 — a PN or a ring nebula around a 
massive WR star? 

The analysis of the spectrum of the nebula of Abell 48 gives 
us indications that this object is most likely a PN, and not 
a nebula around a massive star: 

Diagnostic diagrams are very often used to classify emis- 
sion line sources and separate them from each other on 



the base of the most important line ratios (see e.g., Kni- 
|azev et al.| |2008a| |Frew fc Parker| |2010[ and references 
therein). The most frequently used diagnostic diagram is 
log(F(Ha)/F([Nn] AA6548,6584)) versus log(F(Ha)/F([Sn] 
AA6717,6731)). For Abell 48 these ratios are 0.44 vs. 1.37. 
With these values, Abell 48 is located in the middle of the 
region occupied by PNe and far away from the place where 
Hn regions and WR shells are typically found. 

The density of the Abell 48 nebula and the diagnostic 
diagram logQSll] A6717/A6731)) versus log(F(Ha)/F([Nll] 
AA6548,6584)) also puts this source into the region of PNe 
(cf. Figure[9|, but not of WR shells, which all have densities 



similar to those of H II regions ( Esteban et al 



1992 



Garnett fc Chu||1994| |Greiner et al.|1999||Stock et al 



Fernandez-Martin et al.||2012 i 



1994 



2011 



7.2 Electron-scattering line wings 

Strong emission lines in WR-type spectra show extended 
line wings which can be attributed to line photons that are 
redistributed by scattering on free electrons in the wind. The 
strength of these wings is sensitive to the degree of wind 



inhomogeneities, the so-called clumping (see Section 4.2|. 
The spectra of massive WN stars are typically consistent 
with a clumping density contrast of D = 4 . . . 10 ( Hillier 



1991||Hamann fc Koesterke|19 98 Ham ann et al.|2006 l 

In the spectrum of the CS of Abell 48 these electron- 
scattering line wings are very weak and barely visible, while 
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Figure 8. One-dimensional reduced spectrum of the nebula. Most of the detected strong emission lines of the nebula are marked. All 
detected and measured lines for the nebula are listed in Table [4] The lower spectrum is scaled by 1/30 to show the relative intensities of 
strong lines. 



in contrast the massive WN models display significantly 
stronger wings. 

However, the models which we have calculated for cen- 
tral star parameters, i.e. with logi/L© = 3.7 and M — 
0.6 Mq (cf. Section [5}, show only weak electron-scattering 
line wings with the same values for the clumping contrast as 
for massive WN stars, consistent with the observation. We 
consider this as an evidence for the low-mass nature of the 
CS of Abell48. 



7.3 Location and luminosity of the CS of Abell48 

The best fit to the SED of the CS of Abell48 was achieved 
with a color excess Eb-v = 2.10 ± 0.05 mag (which agrees 
well with Eb-v = 2.15 ± 0.01 mag derived from the line ra- 



massive WN star of logL/L© > 5.3 ( Hamann et al.| 2006), 



tios in the spectrum of the nebula; see Section 6.2 1, implying 
an extinction in the visual of Av ~ 7 mag. The latter trans- 
lates into Ak s ~ 0.8 mag, if one adopts the extinction law 



from Rieke & Lebofsky ( 1985 1. These values should be com- 



pared with the full Galactic V and _?f s -band extinctions in 
the di rection of Abell48, which according to Schlegcl et al. 
( 1998p 1 are ~ 34 and 4 mag, respectively. Thus, one can 
expect that Abell48 is located at a factor of ~ 5 smaller 
distance, d, than the line-of-sight depth of the Galaxy of 
« 20kpc (e.g. |Churchwell et al.|[2009}, i.e. at d w 4kpc. 



Similar distance estimate (d ss 3.9 kpc) could be derived if 
one assumes an average of 1.8 mag of visual extinction per 
kpc in the Galactic plane. 

These distance estimates are supported by a study of 
the distribution of interstellar extinction towards I w 29°, 
which shows that Ak b grows from «0.3 to 1.0 mag between 
2.5 and 4 kpc, then reaches a values of ~1.5 mag at 6 kpc, 
and then abruptly increases to «2.5mag between 6 and 
7 kpc (Negueruela et al. 2011[ ). Since Ak s towards Abell48 is 
~0.8mag, one has that this object should be located in the 



Scutum-Centaurus arm at d < 4 kpc (see fig. 12 in Negueru- 



ela et al. 2011 1 and that its luminosity of logL/Lo ^ 4.3 
is much smaller than the lowest plausible luminosity for a 



but within the range of luminosities derived for CSPNe (e.g 
Miller Bertolami fc Althaus|20"07| ). We adopt the typical lu- 
minosity of CSPNe of L = 6000 Z/q (see Section [5J, which 
leads to a distance to Abell48 of d — 1.9kpc. Note that the 
luminosity of Abell48 and its i?* and M (given in Table |3J> 
can be scaled to different distances as oc d 2 , oc d, and oc d 3 ", 
respectively ( Schmutz et al.|1989 l. 

Alternatively, if the CS of Abell 48 were a massive WN 
star, then it should be located well beyond the Sagittarius 
Arm. Indeed, if using the mean A" s -band absolute magni- 



tude of massive WN5-6 stars from Crowther et al. (20061, 
Mk b = —4.41 mag and the .Kg-band extinction derived from 
the SED fitting and line ratios in the spectrum of the nebula, 
one gets a distance modulus of 15.94 mag, which translates 
into a distance of 15.4 kpc. An even larger distance would be 
derived, if one estimates Ak s from the intrinsic J — K B color 



of WN5-6 stars of 0.18 mag given in Crowther et al. (20061. 
In this case, one gets Ak b = 0.66 mag and d = 16.4 kpc. 
These two distance estimates imply that Abell 48 would be 
located either in the Perseus or the Outer Arm. This loca- 
tion, however, is in contradiction with the moderate extinc- 
tion towards Abell 48 in view of the presence of an extinction 
wall at 



6 kpc (e.g. Negueruela et al.||2011[ ). 



7.4 Abell 48 as a runaway 

The short distance to Abell 48 is also supported by the de- 
tection of a significant proper motion for its CS: \x a c ost? = 
-13.3 ± 4.8masyr _1 and fi$ = — 14.8 ± 4.8masyr _1 ( |Roser| 



et al.||2010[ ). To convert this proper motion and the radial 
velocity of Abell 48 (see Table|3| into the peculiar transverse 
and radial velocities, we used the Galactic con stants Ro — 
8.0 kpc and 9 = 240 km s" 1 ( |Reid et al.||2009[ ) and the so- 
lar peculiar motion (U q ,V q ,Wq) = (11.1, 12.2, 7.3) kms" 1 
(Schonrich et al. 20101. The derived velocity components 



and the total space velocity are given in T able |6| an d im- 
ply that Abell 48 is a runaway system (e.g. Blaauw 19611 



See also http://irsa.ipac.caltech.edu/applications/DUST/ 



For the error calculation, only the errors of the proper mo- 
tion and radial velocity measurements were considered. The 
peculiar radial velocity of Abell 48 contributes only slightly 
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WR nebulae 

G 2.4+1.4 (WR102) 
MR 26 (WR 23) 
S 308 (WR 6) 
NGC 3199 (WR18) 
RCW 104 (WR75) 
L 69.8+1.74 (WR131) 




log(/(Ha)//([SII])) 

Figure 9. Diagnostic diagram Ha/([Sn] AA6717, 6731) inten- 
sity ratio versus [Sn] AA6717/6731 ratio as a measure of electron 
density versus excitation for gaseous nebulae. Shown are the lo- 
cations for H II regions as well as the positions of WR nebulae 
from ||Esteban et al.||1992|l. Most PNe lie on the indicated strip 
from ( |Riesgo-Tirado fc L6pez|2002| . Also labeled are the positions 
of the [WN] and [WN/WC] stars Abell48, IC4663, and PBS to 
demonstrate their CSPNe status. 



4, 
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Figure 10. 2° X 2° MSX 8.3 ^m image of the Galactic plane 
(shown by a dashed line) centered at / = 29°, b = 0?5, with the 
position of Abell48 indicated by a circle. The inset shows the 
IRAC 8 fim image of Abell 48. The arrows on both images show 
the direction of motion of Abell 48. The orientation of the images 
is the same. 



Table 6. Abell 48: Peculiar transverse (in Galactic coordinates), 
radial and total space velocities, and height over the Galactic 
plane for different adopted distances. 



Table 7. Measured line ratios from the spectrum of the CS of 
Abell 48 for spectral subtype classification. 



d 


Vl 


Vb 




^pec 


ds'mb 


[kpc] 


[km s —1 ] 


[km s -1 ] 


[kms" 1 ] 


[kms -1 ] 


[PC] 


1.9 


-149±43 


53±43 


37±4 


162±42 


15 


4.0 


-296±91 


104±91 


-8±4 


314±91 


32 


15.4 


-970±350 


378±350 


74±4 


1043±349 


122 



to its total space velocity, v pcc (Table [6J. These estimates 
show that the larger d the larger v pcc and correspondingly 
the lower the likelihood that the CS of Abell 48 can be accel- 
erated to such high velocity. Consequently, one can conclude 
that the short distance to Abell 48 is more likely. 

The runaway status of Abell 48 is consistent with the 
orientation of its peculiar velocity with respect to the Galac- 
tic plane. Fig.[To| shows the MSX 8.3 /im image of a 2° x 2° 
region of the Galactic plane (centered at I — 29°, b = 0?5). 
One can see that Abell 48 is moving in the "correct" direc- 
tion, i.e. away from the Galactic plane. On the other hand, 
the large margin of error in Vb leaves the possibility that 
Abell 48 is moving almost parallel to the Galactic plane, 
which in turn is consistent with the small height above the 
Galactic plane (cf. Table |6f. 



7.5 Classification 

Based on the arguments above we conclude that Abell 48 is 
a PN with a low-mass central star. 

The abundance pattern, that our spectral analysis re- 
vealed, corresponds to a WN-type composition, where is he- 
lium dominates, nitrogen is enhanced and the carbon and 



line ratio 


peak/ continuum 


N v 4604 /N in 4640 


2.0 


Civ 5808 /He n 5411 


0.7 


Civ 5808 /He 1 5878 


1.7 



oxygen abundances are roughly solar. Therefore, we sug- 
gest to assign Abell 48 to the [WN] class, becoming thus the 
second known member of this spectral class after IC 4663 
( Miszalski et al.|2012 |. 

To determine the detailed subtype of Abell 48, we apply 
the 3D classification scheme by Smith et al. (1996). The 
measured line ratios (see Table |2| and |7j suggest that the 
star belongs to the [WN 5] subtype. 

In this connection, we note that none of the massive 



WNE stars possess a compact circumstellar nebula (Gvara- 



madze et al. 2009| |2010a| and references therein), which 
provides one more evidence in support of the PN status of 
Abell 48. 



7.6 Evolutionary status 

Spectroscopically, Abell 48 with its strong helium and nitro- 
gen emission lines belongs clearly to the [WN] class. How- 
ever, Abell 48 seems to be different in some respect from 



the other [WN] star IC4663. While Misza lski et al.| (2012) 
could not find any residual hydrogen in IC 4663, the pres- 
ence of some hydrogen is indicated by our optical spectrum 
of the CS of Abell 48. Such residual hydrogen would rule 
out a merger scenario with a sequence He- WD + He- WD — > 
[WN] — > O(He) as proposed by Reindl et al. (in prep.) for the 
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O(He) and [WN] stars. Moreover, if the nebula of Abell48 
has been ejected during the merging process it would be 
hydrogen-deficient, contrary to the observation. Thus, there 
must exist different evolutionary channels for the formation 
of these two [WN] stars. 

For the [WN/WC] star PB 



Miller Bertolami et al. 



(20111 suggested a diffusion-induced nova (DIN) of a low- 
mass (M < 0.6 Mq) post-AGB object as the origin of such 
stars. Among the models calculated for the different sce- 
narios of the DIN there is one which includes overshoot- 
ing in the convective zone generated during the CNO-flash, 
which has abundances that are very similar to those that 
we found for Abell48. The surface abundances of the model 
are H:He:C:N:O=22:73:0. 66:3. 8:0.1 per cent by mass, that 
is close to our results of H:He:C:N:O=10:85:0.3:5:0.6 for the 
wind of Abell48. 

The nebula abundance ratio N / O = 0.25 is below 



the sharp limit of N/0> 0.8 that Kaler fe Jacoby ( 1989 1 



deduced for N-enriched PNe, which are supposed to indi- 
cate hot bottom burning in more massive AGB-stars with 
M core > 0.8 Mq. The low core mass of Abell48 would be in 



line with the models by Miller Bertolami et al. (2011 



However, the stellar evolutionary models by Miller 
|Bertolami et al.| ( |201l] ) assume a low-metallicity (Z — 0.001) 
progenitor, which is at least for the CS of PB 8 not the case, 
as we inferred solar-like iron-group abundances from the UV 
spectra of this object. For the CS of Abell48 we could not 
determine any iron-group abundances, but the metallicity of 
the nebula is roug hly solar (cf. Table |Hj|. 

Further problems for this scenario are posed by the in- 
volved timescales. Assuming a PN radius of 0.2 pc as de- 
duced in Section |6.1| and a typical expansion velocity of 
SOkms" 1 , the nebula age is only 6500 yr. This is far too 
short compared to the evolutionary time since the ejection 



of the nebula at the AGB (10 b - 10' yr, Miller Bertolami 



et al.|20lT l, especially since this scenario works only for low 
stellar masses M < 0.6 Mq 
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